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Including collisional decoherence explicitly, phase sensitivity for estimating effective scattering 
strength \ °f a two-component Bose-Einstein condensate is derived analytically. With a measure- 
ment of spin operator J x , we find that the optimal sensitivity depends on initial coherent spin state. 
It degrades by a factor of below super-Heisenberg limit oc l/N 3 ^ 2 for particle number N and 

the dephasing rate 1 « 7 < N 3 ^ A . With a J v measurement, our analytical results confirm that the 
phase (j> — xt ~ can be detected at the limit even in the presence of the dephasing. 



PACS numbers: 03.75.Dg, 03.75.Mn, 03.75.Gg 

Parameter estimation with its precision beyond stan- 
dard quantum limit (SQL) is a long-standing challenge 
in quantum metrology. The achievable precision depends 
on the input state llHZl , the observable being measured 
at output ports Ia4ip|. and the coupling nature of the 
system Hamiltonian [ll| - |l5j . In standard Ramsey inter- 
ferometry, for instance, resonant atomic frequency \ is 
estimated with a coherent spin state (CSS) that evolves 
freely under a linear coupling %J Z . It has been shown 
that the precision scales as 1/iV 1 / 2 (i.e., the SQL) for J x 
or j y measurement [l3| . The phase sensitivity can be 
improved to the so-called Heisenberg limit 1/N with an 
entangled input state Here, collective spin opera- 

tors J v = i <7^ with the Pauli matrices <J v =x,y,z- 

Recently, atom interferometry with Bosc-Einstcin con- 
densates (BEC) becomes a topical area of study due 
to potential applications in quantum metrology [7j and 
quantum information [161 ] . Elastic collision of condensed 
atoms can be described by the 'one-axis twisting' Hamil- 
tonian xJz 12 1; capable of creating multipartite entan- 
glement and spin squeezing [l6l - |20| . Control of spin dy- 
namics requires precise measurement of the effective in- 
teraction strength \ [2l[ . Rey et al. [l2j recently claimed 
that even with an initially CSS, the Heisenberg limit is 
achievable if the free evolution under \ J z is simply re- 
placed by x^ ■ A better scaling oc 1 /TV 3 / 2 is proposed for 
the CSS prepared by a tt/4 or 3tt/4 pulse [II [3. Such 
a super-Heisenberg scaling can also reach in nonlinear 
optical and nano- mechanical systems (Til [l5j . 

In this brief report, we investigate carefully the ef- 
fects of collisional dephasing on the BEC-based quantum 
metrology (T214T4T ]. described by the Hamiltonian (ft = 1): 
H = \J Z + QxJx + QyJy, where the effective interaction 
strength \ an d the Rabi frequencies fl x = Re(0) and 



fly = lm(Q) are tunable in real experiments [19l . |20| |. 
The metrology protocol starts with a Ramsey pulse ap- 



plied to the BEC with all spin up, yielding a product 



CSS: |* 6 



-iOJv 



J, J) = \0,0), where the polar an- 



gle 6 = flyt given by the pulse area. Next, the system 
evolves freely for a time t, \tyg ((/>)) = e~ l ^' J z \^e)i with 
a dimensionless phase shift cf> = xt. Finally, an equato- 
rial component of the total spin J x or J y is measured to 
estimate x 0- 

For nonzero x, the accumulated phase <f> manifests it- 
self as oscillations of the Ramsey signal (J x ) or (J y ), with 
its sensitivity quantified by (l3| 



S(f) v = tSxv 



\d(J v )/d<p\ 



(v = x,or y), 



(1) 



where the variance AA = ((A 2 ) — (A) 2 ) 1 / 2 and the expec- 
tation value (A) = (</>) 1 Xl'e (</>)) ■ The J x / y measure- 
ment is achievable by applying a second 7r/2 pulse after 
the free evolution, e~ l ^ Jy ^ \^g((f>)), and following with a 
detection of population imbalance (i.e., J z ), as done in 
standard Ramsey interferometry. 

Collective spin excitation and external field fluctua- 
tions leads to an enhanced phase diffusion of the BEC 
[22j ] . To describe it qualitatively, we assume that free evo- 
lution of the system obeys the master equation [H, [24| : 
p = i[p, xJ 2 } + T(2J z pJ z - J 2 p - pJ 2 ), where T denotes 
the dephasing rate. For single-particle case, the second 
term reduces to \{a z pa z — p) [25|]. The many-body de- 
coherence considered here, known as the collisional de- 
phasing (23[, can be solved exactly with the density ma- 
trix element p m , n {4>) = (J,m\p\J,n) — p m ,„(0) exp[i(n 2 - 
m 2 )4> — 7(771 — n) 2 4>], which yields 
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where J = N/2, 7 = T/x, and (j) = xt- Note that the 
collisional dephasing imposes an exponential decay to av- 
erage value of J + and its higher moment j?, but main- 
tains that of J z and J 2 intact, i.e., (J z ) = J cosO and 
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FIG. 1: (Color online) Phase sensitivities 5<f> x (left) and 5(f> y 
(right) as a function of phase shift (f> [in units of iy/(^/2J)] for 
J = N/2 = 10 3 and 7 = (up), 10 2 (bottom). Solid (dashed 
red) lines are analytical (exact numerical) results for 9 = n/4 
(a, b, d) and 9 = 7r/6 (c). Empty circles in (a) is exact result 
for 9 — n/2, while in (c) and (d) are obtained from Eq. ([9]) 
and Eq. ()10|) for 9 — tt/6 and 7r/4, respectively. Horizon- 
tal lines denote the Heisenberg limit 1/(V2J) [12], and the 
super-Heisenberg limit l/(y/2J 3 ^ 2 ) Vertical lines in (a) 

and (b) correspond to \<j>\ = 0.89tt/(^/2J) and 0.257r/(v / 2J), 
respectively. 



(J 2 ) = J 2 — J(J - 1/2) sin 2 6> (c.f. Ref. [3). The slope 
of the signal d(J x / y ) /d<f> corresponds to real or imaginary 
part of d(J+) /d(j), with 



d(J 4 



= J(2J- I)e~ 7<?i (icos0cos<?!>-sm0) 

x sin(0)(cos0 + zcos0sin(/>) 2J ~ 2 . (4) 



We exclude the case 9 — and tt, due to d(J v )j/d(j) = 
and thus 8<p v — > 00, which implies no information about 
<f> (or x for a given £) is gained from J x and J y measure- 
ments. 

To obtain scaling rule of 8<p v , we perform standard 
short-time analysis to the above exact results pj], fisl j . 
After some straightforward calculations, we get (A J x ) 2 ~ 
^[1— (170 — ^171) sin 2 0], where rj = ^(l + e -4 ' 3 cos 2a) and 
r/i = (1 - e- 2/3 )(l - e~ 2/3 cos2a) + 2<?!)e~ 2/3 cos(6») sin 2a. 
In the short-time limit (\<f>\ << 1), the parameters a = 
2 J '(/> cos 9 ~ 1 and /3 = J</> 2 sin 2 0+7</> << 1, which enable 
us to expand 770 ) up to the zeroth (the first)-order of 
P. Similarly, Eq. (j4j gives d(j x )/d<f> ~ - J 2 sin (a) sin 20 
for 9 7^ 7r/2. As a result, we obtain 



1 + (cos 9 cot a + 2 Jcj> sin 2 0) 2 + 47 J4> sin 2 
2 J 3 sin 2 20 ' 



(5) 



Replacing a with a + 7r/2, we also get the short-time 
solution of (AJ y ) 2 and that of the sensitivity 



1 + (cos 9 tan a — 2J</>sin 2 



AjJcf) sin 



2 J 3 sin 2 26» 



(6) 



As depicted in Fig. [TJ the sensitivities oscillate rapidly 
in a fringe pattern and diverge at \cf>\ — sir/ (2 J cos9) 
and (s + l/2)7r/(2 J cos 9) for an integer s = 0, 1, etc., 
given by cot a — > 00 and tana — > 00 in Eq. (JS|) and 
Eq. ([6]), respectively. Within central few fringes (s < 2), 
our analytical results (thin solid lines) are coincident with 
the exact numerical simulations (red dashed lines). 

We now analyze the achievable sensitivity for J x and 
J y measurements one by one. Firstly, let us consider the 
case 7 = 0. Via minimizing Eq. ([5]) with respect to 0, we 
find that local minima of 86 r occur when 



cos(0) cot a + 2J<f> sin 2 9 = 0, or sin a 



= cot 9. 

(8<f>x)it 



(7) 



The first transcendental equation gives 
l/CV^./ 3 / 2 ! sin26l|), as predicted in Ref. [U]. For 9 = 
7r/4 or 37t/4, it becomes l/(\/2J 3 / 2 ). Such a super- 
Heisenberg scaling of the sensitivity appears when a = 
y/2J</> ~ ±0.89tt, i.e., <£ min = \x\mint ^ 0.89tt/( v / 2J) [see 
Fig. [31a)], which is valid provided J > 10 2 [Fig. E^a)]. 
Numerically, we also find that the best sensitivity phase 
</>min depends on polar angle of the initial state 9 [empty 
circles of Fig.^b)]. For 7r/4 < 9 < 37r/4, it is in fact de- 
termined by sin a = cot 9, which predicts m in — 0.2tt/J 
for 9 = 7r/3 [crosses of Fig. HJa)]. 
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FIG. 2: (Color online) The best-sensitivity phase m in as a 
function of J (a) and 9 (b) for measurement and 7 = 0. 
In (a), numerical results of </> m in for 8 = 7r/4 (empty circles) 
and 7r/3 (crosses), fit with 0.897r/(v'2J r ) (dashed red line) and 
0.2n/J (solid line). In (b), contour plots of the first (dashed 
green line) and the second (solid blue line) equations of © 
for J = 10 3 , with numerical m j n (empty circles). The shade 
region denotes the central fringe < (frmin < tt / (2 J cos 6). 



To proceed, we consider the J y measurement in the 
absence of the dephasing. From Eq. ©, one can 
find that minimal value of the sensitivity (<50 y ) m in — 
1/(V2J 3/2 | sin20|) occurs at ~ 0. Obviously, the 
super-Heisenberg limit 1/(a/2J 3 / 2 ) is attainable for the 
optimal CSS with = tt/4 or 3tt/4 [see Fig. [TJb), also 
Ref. [3). 

The above results can be casted in a more transparent 
form by setting (S(f> v ) m i n = nJ~^ v (for v = x, y) [l3j . 
with a pre-factor k and scaling exponent 



& = - 



ln[V2(^) r 



ln(|sin20|) 



(8) 



In J 2 In J 

where we set k = 1/V2 to ensure — > 3/2 as — >• 
7r/4 or 37r/4. One can find from Fig. [3] that without the 
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FIG. 3: (Color online) Scaling exponents £ x and £ y as a func- 
tion of polar angle 6 for J = 10 5 . The circles (the crosses) 
are numerical results of £ x (£ y ) for 7 = 0, and the dashed 
red line is obtained from Eq. @. The triangles are numerical 
results of £ x for 7 — 10 3 , and the solid line is predicted by 
Eq. (|lip . Vertical lines denote 9 = 7r/6, 7r/4, 37r/4, and 57r/6; 
horizontal ones correspond to 3/2, 3/2 — ln(27)/(3 In J). 



dephasing, numerical results of £ x (empty circles) and £ y 
(green crosses) agree with Eq. ([5]) (dashed red line). 

Finally, we discuss the achievable sensitivity in the 
presence of the dephasing. For 9 ^ ir/2 and 7 >> 1, 
Eq. ([5]) and Eq. ([6]) can be simplified as 



1 



[(4J(/>)- 2 +4jJcj) sin 2 6] , (9) 



2 J 3 sin 2 26> 

^ - ijd^^^- 2 ^], do) 

which correspond to the envelope curves of the sen- 
sitivities [see empty circles of Figs. [TJc) and (d)]. 



For J x measurement, the best sensitivity 
37 2 / 3 /(8J 3 sin 2/3 6>cos 2 9), and thus 



lnC^^- 1 / 3 sin 1/3 (9 Icos^l) 



x /min 



fa - 2 



In J 



(11) 



As shown in Fig. [3l our analytical result shows a good 
agreement with numerical simulations (triangles) for the 
dephasing rate 1 « 7 < J 3 / 4 . Remarkably, the optimal 
scaling t; x ~ 3/2 — 111(27)/ (3 In J) i s obtained at 9 = 
it/ 6 or 57r/6, which leads to the best-sensitivity phase 
min ~ l/[J(2 7 ) 1 /3] W ith (5<fe) min ~ (2 7 ) 1 / 3 /(72J 3 / 2 ). 
For J y measurement, the phase |0| ~ sn/(y2J) (with an 
integer s < 2) can be detected at the super-Heisenberg 
limit provided 9 — ir/4 or 37r/4. However, the sensitivity 
for |0| > 27r/(v / 2J) degrades rapidly [see Fig.QJd)]. 



In summary, we have derived analytical results for pre- 
cise estimation of effective interaction strength \ in a 
two-component BEC. Without collisional dephasing, the 
best-sensitivity phase (f> m - m for J x measurement depends 
on the initial CSS \9, 0) and bifurcates at 9 ~ 7r/4 or 37r/4 
(see Fig. [5}d). In the presence of the dephasing (with 
1 < 7 < J 3 / 4 ), the optimal CSS becomes \9 = tt/6,0), 
and the achievable sensitivity is reduced by a factor of 
(27) 1 / 3 below the scaling l/(\/2J 3/2 ). Our analytical 
results confirm that the detection of </, 

surement to the output state e "2-*e T " ! e 
shows its advantages since phase estimation of 4> ~ can 
reach the super-Heisenberg limit even in the presence of 
the dephasing [l3j . 



y , i.e., the J z mea- 
e-^e-^^lJ, J), 
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